RATIONAL MATRIX SOLUTIONS TO THE LEECH EQUATION: 
THE BALL-TRENT APPROACH REVISITED 



SANNE TER HORST 



Abstract. Using spectral factorization techniques, a method is given by which 
rational matrix solutions to the Leech equation with rational matrix data can 
be computed explicitly. This method is based on an approach by J. A. Ball 
and T.T. Trent, and generalizes techniques from recent work of T.T. Trent for 
the case of polynomial matrix data. 



0. Introduction 

Consider ff°°-matrix functions G G H™ Xp and K G H™ Xq , and let T G : £ 2 + (CP) ->■ 
£\(C m ) and T K : ^(C 9 ) t\(C m ) be the corresponding (block) Toeplitz oper- 
ators. See Section [T] below for the definitions of these spaces and operators. A 
beautiful unpublished result of R.B. Leech (cf., [11]) tells us that there exists an 
X G H^ n such that 

(0.1) G(z)X{z) = K{z) (zeB), and ||JC||oq < 1, 

with P the open unit disc in C, if and only if 

(0.2) T G T G - T K T* K is positive. 

Note that (jO.ip is equivalent to TqTx — Tk and ||Tx|| < 1. Hence Leech's theorem 
can be viewed as the analogue of the Douglas factorization lemma [7] within the 
class of analytic Toeplitz operators. The necessity of (|0.2|) follows directly from 
Douglas' factorization lemma and the reformulation of ([O.ljl in terms of Toeplitz 
operators. The other implication is more involved. The solution criterion (|0.2[) can 
also be formulated directly in terms of the functions G and K, it is equivalent to 
the map 

, no . Tt s G(z)G(w)* - K(z)K(wY . 

(0.3) L(z,w) = , = LJ — (z,weB) 

1 — zw 

being a positive kernel in the sense of Aronszajn [I], that is, for any finite sequence 
Zii • • • ) z n G ID> the block operator matrix [L(z^ 2j)]ij=i,..., n defines a positive op- 
erator on the Hilbert space direct sum of n copies of C m . We note that the actual 
result by Leech is stated in the general context of Hilbert space operators inter- 
twining shift operators, and in particular holds for operator-valued iJ°°-functions 
as well. Our interest is primarily in the case where G and K are rational matrix 
functions. 

There exists various proofs of Leech's theorem, see [10] and the references therein. 
In [3] Ball and Trent prove a generalization of Leech's theorem to the polydisc in 
C d , adapting a technique coined the 'lurking isometry' approach in [2], and give a 
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description of all X E H^ q satisfying ()0.1[) . We briefly outline the construction 
here, specified to the single variable case. 

The positivity of T G T£ - T K T^ implies we can factor TqTq - T K T£ as TqT g - 
T K T^ = A A*, for some operator A Q : H -> £\{C m ) such that Ker A = {0}. The 
latter implies that dim% = rank [TqTq — TrT^). Such a factorization is often 
referred to as a Kolmogorov decomposition in the literature, cf., [6]. Let A be 
the analytic operator- valued function on D, with values A (z) : H — > C m , zeD, 
defined by A (z)h = ($ m A h)(z), zeD,/i£ % . Here g" m is the Fourier transform 
mapping ^(C m ) isometrically onto the Hardy space H^. Next one verifies that G, 
K and A satisfy the following identity: 

zwk {z)k {w)* + G(z)G(w)* = 

(0.4) = A (z)A (w)* + K(z)K(w)* (z,w G B). 

From this identity one derives the existence of a partial isometry 

^ M °=[t D 

such that 

(0.6) [ zA (z) G(z) )M =[ A (z) K(z) ] (z G D). 

This in turn implies that the function X defined on B by 
(0.7) X(z) = D + zC (I-zA )- 1 B (z G B) 

is in H^. and satisfies (|0.1[) . If one considers all contractions M of the form (|0.5j) 





7^o 




Ho 






— > 





such that (|0.6p holds, possibly enlarging Ho, all solutions X to (|0.ip are obtained 
via (|0Jl) . 

From the point of view of rational matrix functions the above construction has 
one disadvantage. In general, the Hilbert space T-L appearing in (|0.5|) is infinite 
dimensional, and in that case it is hard to see when the solution X in (|0.7[) is 
rational. In fact, even if both G and K are rational matrix functions, TqTq — TkT^ 
may very well be of infinite rank. More precisely, see Theorem 13.21 below, in the 
rational matrix case T g Tq — T K T^ has finite rank if and only if G(e u )G(e lt )* = 
K(e lt )K(e it )* for all < t < 2tt. Overcoming this difficulty is the main theme of 
the present paper. 

In the context of the Toeplitz-corona problem, which can be reduced to the 
special case of (jO.ll) with q — m and K(z) — I m , z G B, Trent [12] deduced a 
modification of the above procedure for the special case that G is a row vector 
(m = 1) polynomial, leading to a rational column vector solution of McMillan 
degree at most the highest degree of the polynomials occurring in G. Throughout 
this paper, the McMillan degree of a rational matrix function V will be denoted 
by S(V); see Section Q] for the precise definition of 8(V). The procedure of [12] was 
recently extended in [13] to the general case of the Leech equation (|0.1[) . with G 
and K rational matrix functions, by reducing it to the case where G and K have 
polynomial entries, and solving the latter problem via techniques similar to those 
in [H]. 

In the present paper we also consider the Leech equation (jO.ip with G and K 
rational matrix functions. However, instead of reducing to the case of polynomial 
data, we associate our problem with another Leech equation, with data functions G 
and K, i.e., with the same G. The advantage of our approach is that we keep better 
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track of the McMillan degrees in our computations, leading to sharper bounds on 
the McMillan degrees of the solutions. The construction of K even works in the 
case where G and K are not rational, provided that the function R £ L™ xm defined 

by 

(0.8) R{e lt ) = G(e lt )G(e lt )* ~ K(e u )K(e u )* (a.e. t £ [0, 2tt]) 

admits an outer spectral factor, that is, a function $ £ H^ xm , for some r < m, 
with Tr = T£Tg> and kerT^ = {0}. Note that outer spectral factors are unique 
up to multiplication with a unitary constant matrix on the left, hence, with some 
abuse of terminology, we will refer to the outer spectral factor, provided it exist. If 
G and K are rational, then so is R, and this implies an outer spectral factor of R 
exists. 

Our method requires the following procedure: 

1. Define R £ i™ xm by (|0.8|) . Then Tr is positive. Assume R admits an 
outer spectral factor $ £ H^ m , for some r < m. 

2. The subspace 

(0.9) Mi := {/ £ 4(C r ) | T|/ £ ImH G + lmH K }. 

is invariant under the backward shift on £^_(C r ), and hence, by the Beurling- 
Lax theorem, there exists an inner function £ H?t k , for some k < r, such 
that the range of Tq is the orthogonal complement of A4<s>. 

3. Define F £ L™ xk by F{e lt ) = $(e lt )*6(e"), for a.e. t £ [0,2tt]. Then 

The claims in the above steps will be proved Section [5] The function F defined in 
Step 3 can be taken as a particular choice for the function F appearing in the next 
theorem. This theorem provides the basis for our method and is the main result of 
the present paper; a proof will be given in Section [2j 

Theorem 0.1. Assume G £ H™ xp and K £ H£ xq such that T G T£ - T K T* K is 
positive and the function R defined in (|0.8|) admits an outer spectral factor. Then 
there exists a function F £ H^ xk , for some k < m, such that: 

(i) TgT'q — TkTk — TpTp is positive; 

(ii) rank (T G T£ - T K T* K - T F T F ) < dim (lmH G + ImH K ^ . 

Here Hq and Hk denote the Hankel operators of G and K , respectively. 

Given F as in Theorem I0.1[ we apply the Ball- Trent approach with K replaced 
by K = [K F]. This yields tf^-solutions X = [X Y] of 

(0.10) G(z) [X(z) Y(z)] = [K(z) F(z)] (\z\ < 1), and || [X Y] |U < 1. 

Note that (|0.10|) implies that X satisfies (|0.1|) . Whether or not all solutions of (|0.1j) 
can be obtained via this procedure is still an open problem. 

This procedure is specifically of interest in case G and K are rational matrix 
functions. In that case the upper bound in (ii) is finite, and serves as an upper 
bound on the least possible McMillan degree of solutions X to (|0.10j) . hence the 
same upper bound applies to X. The following theorem provides some additional 
results for the case of rational data functions; a proof will be given in Section |3] 

Theorem 0.2. Let G £ H^ xp and K £ H^ xq be rational matrix functions such 
that TqTq—TkFk is positive. Then the function R £ i™ xm defined by (|0.8j) admits 
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an outer spectral factor $. Moreover, in this case the functions R, $, and F 
defined in the above procedure are all rational matrix functions whose McMillan 
degrees satisfy 

(0.11) -8(R) = ($(*) < 5(0) = 6(F) = dimM$ < oo, 

and 9 is two-sided inner, i.e., k = r and Q(e lt )*Q(e lt ) = I r = 6(e I ')6(e 1 ')* for 
each t £ [0, 2tt]. Finally, we have 

fCi 1 e)\ *F *F* *F *F* *F *F* / / 7T* i / / TT* 1 1 TJ* 

(U.iZ) IqIq — 1k±k~ 1 F^f — tl Ktl K + tlptlp — lio-tlc- 

In particular, the left hand side in inequality (ii) in Theorem ] U.l\ is equal to rank (HkH^+ 

/ / 7T* / I TT * \ 

liptlp — UqUq). 

Thus, in case G and K are rational matrix functions, the problem reduces to 
computing a Kolmogorov decomposition of the right hand side of f)0. 12[1 . Note 
that there are effective ways to computing Kolmogorov decompositions, cf., [6]. 
Moreover, the functions R, $, G and F can be computed explicitly using state space 
techniques from mathematical systems theory (cf., [H [5]), starting from a state 
space representation of the function [ G K ) . This will be the topic of a forthcoming 
paper of the present author together with A.E. Frazho and M.A. Kaashoek. 

The paper consists of 0] sections, not counting the present introduction. Section 
[T] contains some of the notations and terminology as well as some operator theory 
preliminaries used in the sequel. The main result, Theorem lO.il is proved in Section 
[2] In Section [3] the focus lays on the case that G and K are rational matrix 
functions; a proof of Theorem 10.21 will be given as well as a criterion for the case 
that TqTq — TkT^ has finite rank. The final section contains some general operator 
theoretical results, and their proofs, that are used in the preceding sections. 

1. Preliminaries 

In this section we introduce notations and terminology used throughout the 
paper and we present some operator theory preliminaries. 

With operator we mean a continuous linear map acting between two Hilbcrt 
spaces. In particular, all operators in this paper are by definition bounded. In- 
vertibility of an operator means the operator has a bounded inverse. Let H be a 
Hilbert space. A subspace of H is a closed linear manifold within H. The identity 
operator on T-L will denoted by I-u and the k x k identity matrix by Ik ■ Often these 
subscripts R and k will be omitted. We say that an operator T on R. is positive 
whenever the inner product (Tu, u) > for each u £ H, and T is said to be positive 
definite whenever T is both positive and invertible. The notations T > and T > 
will be used to indicate the positivity, respectively positive definiteness, of T . In 
case T\ and T2 are selfadjoint operators on R, we will write T± > T2, resp. T± > T%, 
to indicate T 1 ~T 2 > 0, resp. T 1 ~T 2 > 0. 

The symbol H^ xp will indicate the Hardy space of all uniformly bounded ana- 
lytic m x p matrix-valued functions in the open unit disc. For any V G ff^ xp the 
supremum norm of V is defined by ||^||oo = su P| z |<i ll^(- z )l| 5 making H^ xp into 
a Banach space. Here we follow the convention that the norm \\M\\ of an m x p 
matrix M is equal to the norm of the operator from C p into C m induced by M in 
the canonical way. We write I^ xp for the Banach space consisting of all Lebesgue 
measurable, essentially bounded m x p-matrix functions on the unit circle T to- 
gether with the essential supremum norm, also denoted by || ||oo. The space H^ x 
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will be viewed both as a sub-Banach space of L™ xp and as a Banach space in its 
own right. 



With a function Z G L™ xp we associate the functions Z' 



defined by 
(1.1) Z*(e u ) = Z{e lt Y 



e L™ xm and Z* G i^ xp 



and Z l {e lt ) = Z(e~ u ) (a.e. t G [0, 2tt]) 



For y G H^ x „, the functions 1^* and V* can be uniquely extended to bounded 

, infinity included, via the formulas 



analytic functions on the open exterior disc 
V*{z) = V(l/z)* and V*(z) = V{l/z), \z\ > 1. 

By £ 2 (C k ) and £ 2 _(C k ) we denote the Hilbert spaces consisting of bilateral, 
respectively unilateral, square summable sequences with values in C . Viewing 
£a_(C) as a sub-Hilbert space of £ 2 (C k ), we write £ 2 _(C k ) for the orthogonal com- 
plement of £ 2 + (<C k ) in £ 2 (C k ). The symbol Sk stands for the (block) forward shift 
on £ 2 + (<C k ), and E/. denotes the canonical embedding of C fc into £+(C k ) defined by 

E k u = [u • • •] T . Note that I - S k S* k = E k E* k . 

Let Z be a function in L^ xp and denote the Fourier coefficients of Z by ... , Z—i, Zq, Z\, Z2, 
Then we define the (block) Toeplitz operator Tz and (block) Hankel operators Hz,+ 



and Hz.- associated with Z by the operators mapping £^_(C P ) into t\{C r ' 
by their infinite block matrix representations 



given 



T z = 



Zq Z-i 

Z-2 Z x 



Z- 

z_ 

Zo 



Hz.+ — 



Z\ Zi z^ 
Z<2 Z-$ z^ 
Z$ Z4 Z5 



Hz- = 



Z- 

z_ 
z_ 



Z- 

z_ 
z_ 



We shall refer to Hz,+ and Hz,- as the analytic, respectively anti- analytic, Hankel 
operator associated with Z. Note that Tz* = T% and H* z , = Hz*-- For V G 
Hmxp we have -f/y,_ = 0, and we will simply write Hy for Hy,+- 

Now consider [/ G H^ xp , V G H™ xp and W G H^ xg . Then the following useful 
identities apply (cf., [5] Proposition 2.14]): 



(1.2) 



Tyw = TyTw, 
Hyw,+ = TyHy/ 



±UV — 1 U 1 V + "17-"V> 
1 -nt/V*,+ — tlijlyt. 



The sets of rational matrix L^< m - and 7?^ m -functions will be denoted by 
9^mxp an< i SK^mxp, respectively. For a m x p rational matrix function Z the 
McMillan degree is denoted by <5(Z) and equals the sum of the local degrees, 
5(Z) = ^2 we cS(Z,w). Here the local degree 8{Z,w) of Z at u> is defined to be 
the rank of the Hankel operator defined by the negative Fourier coefficients of the 
Fourier expansion if Z in a deleted neighborhood of w. See Section 8.4 in [1] for more 
details. It is well known that for V G D\H^ xp the MacMillan degree S(V) equals 
the rank of the Hankel operator Hy 
degree S(Z) is equal to rank (if ^ )+ ) + ra,uk(Hz,-) 



Moreover, for Z G SHL~ XJ) the MacMillan 



2. Proof of Theorem 10.11 

Let G G H^ xp and K G #~ X(? , and define i? G L~ xm by (|01i|) . Throughout 
this section we shall assume that TqTq > TrT^. This implies that R is positive 
on T. Indeed, note that the positivity of the kernel L in (|0.3j) implies that (1 — 
|z| 2 )L(z, z) = G(z)G(z)* — K(z)K(z)* is positive for each z G D. Hence the same is 
true for the non-tangential limits of (1 — |2;| 2 )L(z, z) to the unit circle, which exist 
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for almost all points on the unit circle, where the values coincide with the values of 
R. 

Since the function R is positive on T, it follows that Tr is a positive operator 
on £5_(C m ). Under some additional constraints on the positivity of Tr implies 
that R admits an outer spectral factor (see [9l Proposition V.4.2]), that is, there 
exists a function $ £ H!^ m , for some integer r < m, such that 

(2.1) i? = $*$, i.e., T fl = T|T # , and KerT| = {0}. 

The latter condition says that T$ has dense range, i.e., $ is outer. The function 
$ is unique up to a unitary constant matrix on the left, that is, if \1/ is another 
outer function satisfying R = $* 1 J r , then $ and $ are matrix functions of the same 
size, and $(•) = If^(-) where U is a constant unitary matrix. With some abuse of 
terminology, we shall refer to $ as the outer spectral factor of R. See [11] for 
further details. 

We start with a few preliminary results. 

Lemma 2.1. Let $ be the r x m outer spectral factor of the function R given by 
(10.8p . Set A/$ = ImiJc; + Imi7/f, and let A4$ 6e i/ie inverse image o/A/$ under 
the map T|, i.e., 

(2.2) X* = (T!)- 1 [7V$] = {/ e e 2 +(C r ) I T|J e ImiT G + ImffK}. 

Then A4$ is a subspace o/£?_(C r ), dim.M$ < dimA/$, and .M$ is invariant under 
the backward shift S*. Moreover, 

(2.3) Im H$E m = Im S*T$E m £ M®. 

Proof. Since T| is a continuous linear map, the inverse image of the closed linear 
manifold W$ under T| is again linear and closed. Thus is a subspace. The 
bound on dim Ai<s> follows from the injectivity of TJ. The fact that S^Hg = HgS p 
and S* u Hk = HkSq implies that 

(lmH G +lmH K ) C S* n lmH G + S^lmH K 

= lmH G Sp + Im H K S q C ImTTc + ImiJ A -. 

Thus A/$ is invariant under S^. Take / € A^$, i.e., T£f £ 7V$. Using .§ r T$ = 
T$S m we have 

T£S*f = S^T^f £ C A/*. 

Thus £*/ G Hence .M$ is invariant under the backward shift S*. 

Next we prove (|2.3I) . Inspecting the first columns in and T$ yields H$E m = 
S*T$E m . Hence the identity in ((23} holds. Take u e C m , and put x = S*T®E m u. 
Then 

T$x = T^S*T^>E m u = S* n T^T^E m u — S^TfiEmU 

= ^{TgTq + H G H G )E. m u ~ S^TkT^ + H K H^)E m u. 

Note that T G E m = E p G(0)*, S^T G E p = H G E p and S* m H G = H G S p . Hence 

S* m (T G T G + H G H G )E m u = H G (E p G(0)* + S p H G E)u £ Im H G . 

Similarly, S* n {T K T^ + H K H* K )E m u £ ImH K . This shows that Tgx belongs to 
ImH G + ImH K C A/$, and thus x £ M$. Hence lmS*T 9 E m C M$. □ 

Corollary 2.2. Let $ be the r x m outer spectral factor of the function R given by 
dQUJ. Define by J2J2]). Then ImH 9 C M$. 
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Proof. By (|2.3|) . we see that the range of the first block column of Hq> is in Ad®. 
Since H^S m = S*H<j>, it follows that H$S l m = S* l H<^ holds for any positive integer 
I. The fact that M$ is invariant under S* then shows that for any positive integer 
I 

lmH$S l m E m = lmS* r l H^E m = S; l lmH$E m C S?M* C M*. 

This shows that the range of each column of is in A^$, and thus the range of 
-ff$ is included in A^$. □ 

By the Bcurling-Lax-Halmos theorem, the fact that the space A4<p is invariant 
under the backward shift implies M§ = KerTg for some inner function 9 6 H^ k , 
with k some nonnegative integer, k < r. This O is unique up to a constant unitary 
matrix from the right. Despite this mild form of non-uniqueness, we shall refer to 
© as the inner function associated with the space M®. 

Proposition 2.3. Let $ be the r x m outer spectral factor of the function R given 
by (|0.8[) . and let be the r x k inner function associated with the space M.® in 
(|2.2|) . Then F — $*0 belongs to H^ yk . Moreover, we have 

[}) 1 G 1 G ^ 1 K 1 K ^ 1 F 1 F — — tl G tl G ' U KJ^k — U > 

(ii) rank (T G T G - T K T£ - T F T* ) < dim(I mj ff G + ETHk). 

-ffere Pm* * s the orthogonal projection on £ 2 i _(C k ) with range A4$. If in addition 
H G H* G - H K H* K > 0, then rank (T G T G - T K T* K - T F T*) < dim M<s>. 

Proof. Since $* and are matrix- valued I/°°-functions, we have F G L™ xfe . To 
see that F £ H^xk it suffices to show that Hp,- = 0. However, this is the same as 
showing that H F *,+ = 0. Note that kerT^ = i 2 + {<C r ) Q M g ,k, by definition of 0. 
Hence H&Tq = 0, by Corollary l2~2l Thus the third identity in (fTT2]) yields 

Hf\+ = Pe**,+ = T^H<i> = 0, 

and it follows that F 6 H^xk' as claimed. 

Next we deal with item (i). Since Aiq, = KerTQ and is inner, M.^ = ImTe 
and TqTq is the orthogonal projection onto M$ . In particular, / — Pm® — TqTq. 
Applying the second identity in (11.2[) yields 

(2.4) T R = (T G T G - T K T* K ) + (H G H G - H K H* K ). 

With (HU) and I - P M<S , = T e T£ we obtain 

(T G T G - T K T* K ) + (H G H* G - H K H* K ) =T R = T|T$ = 

Tr# 7 > rr\ I rri* ( j I ) \ rr\ rri* I ) rri , rri* rri rri* rri 

$ ".M* J * + 1 F-L F - 

Here we used that T|Te = T$*Te = Tp. This proves the identity in (i). 

To show T^Pm^Tq—HqIIq+HkHk lii positive and to prove the rank constraint 
on this operator, we apply Lemma 14.11 with the following choices of spaces and 
operators: 

V = 4(C m ), Vi = 7V$ , V 2 -V6^ X = H G H* G - H K H* K , 
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Here A/$ and M$ are the spaces denned in Lemma 12. II In particular, 

XV = Im (H G H G - H K H* K ) C (lmH G + lmH K ^ = jV$ = Vi, 
F-^Vi] = (Tl)- 1 ^] = M$ = Wi. 

Furthermore, we have 

YY* - X = T*T^ - H G H* G + H K H* K = T R - H G H* G + H K H* K 
= T G T G - T K T K > 0. 

Hence T^P M ^T$ - H G H G + H K H* K = YP Wl Y* - X is positive by gZQ), and the 
rank constraint (ii) follows from Lemma 14.11 as well. 

Moreover, note that H G H G — HkHk > translates to X > 0. Thus, by the 
last statement of Lemma 14.11 we find that 

rank (T*P M ^ - H G H* G + H K H* K ) = rank (YP Wl Y* - X) < dim Wi, 

which, together with dimWi = dim.M$, proves the last claim. □ 

We will now prove the main result of the present paper. 

Proof of Theorem 10.11 Let M.$> and A/$ be as in Lemma [2.11 and define F as 
in Proposition 12.31 Thus F = $*9, where is the r x m outer spectral factor of 
the function R, and 6 is the r x k inner function associated with the space A4$ 
in (E2J). We know that F E H% yk . With this choice of F, Proposition O tells us 
directly that items (i) and (ii) in Theorem 10.11 are fulfilled. □ 

3. The case where G and K are rational matrix functions. 

Let G E fRH£ xp and K E fRH™ xq such that T G T* - T K T K > 0. The aim of 
this section is to prove Theorem 10.21 In addition we will derive a criterion for the 
case that rank (T G T G - T K T* K ) < oo. 

In the previous section we observed that Tr > 0, where R E L" xm is given by 
(|0.8|) . Since G and K are rational, so is R, and this, together with Tr > 0, implies 
R admits an outer spectral factor $ E D\H^ rn , sor some r < m, see [TTJ Section 
6.6]. Also note that 5(G) = rankiic < °o and 5(K) — r&nkFtK < oo imply that 
the subspace A/$ of Lemma f2.ll is finite dimensional, and hence the subspace Ai<s> 
in (|2.2[) is finite dimensional, since kerT| = {0}. Then Theorem 4.3.2 in [8] yields 
that the inner function 6 associated with M.<$ is a two-sided inner rational matrix 
function, that is, 6 E 9\H^ r and 69* = 6*9 is identically equal to I r . 

The next proposition provides the relations between the McMillan degrees given 
in Theorem 10.21 

Proposition 3.1. Let G E V\H% xp and K E 9iH™ Xq with T G T G - T K T^ > 0. 
Define as in ()2.2|1 . Then the functions R, $, 6 and F defined in Section^ are 
all rational matrix functions and the following bounds on their McMillan degrees 
apply: 

(3.1) ^S(R) = < 5(F) = (5(6) = dim-M*. 

Moreover, we have HqHq = Pm® an d HpHp = T^Pm^Tq . 
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Proof. Corollary 12.21 implies Imff$ C Hence 

<5($) = rank(iJ$) = dimlmil$ < dim.M$. 

Moreover, we have -£Tr,+ = = T^Hq, by the third identity in (ll.2j) applied 

to V*W = Since $ is outer, KerTJ = {0}, and therefore rank7J fl . + = 

rankiJ$ = 5($). By Tr > 0, we have Hr,- = H R+ . In particular, rankiJR._ = 
rankiTjj , = rankH Rt+ , and thus 6(R) — 2rankiJ^. + = 2(5(<£). 

The fact that is inner with — Ker Tg implies TeTg = I - Pm#- Since 
is two-sided inner, we have 00* = 0*0 = I r , hence Tee* = I- Now apply the 
second identity of (|1.2j) . This yields 

II' tT* 7^ <~p* j rp rp* [ } 

iiQtiQ — Jee* — -Le-Le — 1 ^ J o J e — *\M#< 

Hence 

(5(0) = ranki?e = rank (HqHq) = rankPx^, = dim 
Recall that T = $*0. Hence, by the third identity of (|1.2[) . we obtain that 
Hp = iJ$»e = T|i?e- Since $ is outer, we have KerTJ = {0}, which implies 
5(F) = rankH F = iank(T*H e ) = ranki/ e = 6{Q). Finally, H F = T*H e together 
with H e H^ = P M ^ implies H F H* F = T*P M ^. □ 

Note that dim(Imi? G + lmH K ) < S(G) + 5(K) < oo, since G G 9tff~ xp and 
K E 9\H^ xq . Hence, replacing K by K = [K F], reduces the original Leech 
equation (jO.ll) to one where 

(3.2) rank (T G T* G ~ T K T* K ) < oo. 

We will next focus on the case of the Leech equation where the rank constraint 
(|3.2|) holds. The following theorem provides necessary and sufficient conditions for 
(|3~2|) to hold. 

Theorem 3.2. Let G 6 9iff~ Xp and K G ^H^ xq with T G T G -T K T* K > 0. Define 
R G 9tL^ xm &y f|0 - 8f) . TTien i/ie following statements are equivalent: 

(i) rank (T G T G - T K T£) < oo; 

(ii) T R = 0; 

(iii) G(e B )G(e tt )* = K(e lt )K(e lt )* (t G [0, 2tt]). 
Moreover, in this case 

(3.3) TqT g - T K T* K = H K H* K - i? G JI G 
and 

(3.4) (5(G) < 6(K), 5(K) - 5(G) < rank (T G T G - T K T* K ) < 6{K). 

Here 5(G) and 5(K) denote the McMillan degrees of G and K, respectively. 

Proof. Note that (iii) is equivalent to R(e rt ) — for each t 6 [0,27r], hence to 
T R = 0, since R G «T^ xm . Thus (ii) & (iii). 

The fact that G and K are rational matrix i7°°-functions implies that Hq and 
Hk have finite rank, and thus r&nk(H G H G — HkH^) < oo. From formula (I2.4[) 
it then follows that (i) holds if and only if rankT^ < oo. However, R is a rational 
matrix function with no poles of the circle, and thus continuous on the circle. This 
implies that rankTf?, < oo holds if and only if R(e lt ) — for all t G [0, 2tt], and thus 
T R = 0. Hence (i) (ii). 

The combination of Tr = and formula (|2.4|) gives (|3.3j) . 
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Note that for any positive Hilbert space operators Z and Y on V, the inequality 
Z > Y implies rankZ > rankY. Indeed, by Douglas' Factorization Lemma there 
exists a contraction Q on V such that Y? = QZ? . Hence 

rankY = rankY^ — ra,nk(QZ?) < rank(Z^) = rank(Z). 

Applying this inequality with Z — FLkFL* k and Y = H G H G and noting that Z — Y = 
H K H* K - H G H* G = T G T G - T K T^ > 0, we obtain 

5{K) = rank {H K H* K ) = rank (Z) > rank (Y) = rank (H G H G ) = 5(G). 

If we take Z = H K H* K and Y = E K H* K - H G H* G = T G T G - T K T* K , then clearly 
Z > Y, and thus 

8{K) = rank {H K H* K ) = rank (Z) > rank (Y) = rank {T G T G - T K T* K ). 
In addition to Z and Y, set V = H G H G . Then Z = Y + V implies 

rank (Z) = rank (Y + V) < rank (Y) + rank (V). 
Since 6(G) = rank (V) and S(K) = rank (Z), the last part of holds. □ 

Remark 3.3. If the matrix iJ°°-functions G and K are continuous, then the first 
part of Thcorem l3.2| goes through in a slightly altered form. One only has to replace 
(i) by: T G T G — TkTk ^ s compact. The argumentation is similar to the one given in 
the proof of Theorem 13. 2\ where we now use that H G and Hk are compact, since 
G and K are continuous, and that R being continuous together with Tr compact 
implies Tr — 0, and hence R = 0. 

How restrictive condition p.2[) can be becomes evident when considering the 
Toeplitz corona problem. 

Corollary 3.4. Let G £ ViH™ xp such that T G T G > I, i.e., {(OJ holds with K(z) = 
Im for each z G ID). Then rank (T G T G — I) < oo holds if and only if G is a constant 
matrix function whose value is a co-isometry. 

Proof. Clearly if G is a constant matrix function whose value is a co-isometry, 
then T G T G — I, and hence T G T G - I has finite rank. 

Conversely, assume T G T G — I has finite rank. By Theorem 13. 2\ R = GG* — 
I m = 0. Thus GG* — I m . In particular, the values of G are co-isometries. By 
the second identity in {T2J we have Ip (£m) = T GG » = T G T G + H G H G . Thus 
—H G H* G = T G T G — I > 0. This can only occur if H G — 0, i.e., if G is constant 
matrix function. □ 

Corollary 3.5. Let G E <K#~ xp and K e 9UI™ Xq withT G T G -T K T^ > 0. Define 
R, $ ; 0, and F as in Section^ Then 

(3.5) $*$ = R = FF* and $ = OF*. 

Moreover, Tr > if and only $ is invertible outer, that is, r — m and $ has 
an inverse in FF^ xm . In this case F is invertible in LJJ xm with an anti-analytic 
inverse. 

The first two identities in (|3.5|) say that $ is a right and F a left spectral factors of 
R. The last identity, together with two-sided inner, provides a Douglas-Shapiro- 
Shields factorization of $, cf., Chapter 4]. 
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Proof of Corollary 13.51 The identity <J>*$ = R holds by definition of $. Apply- 
ing Theorem 13.21 with K replaced by K = [K F], where we note that condition (i) 
is satisfied by Theorem lO.il yields 

GG* = KK* = KK* + FF* , i.e. FF* = GG* — KK* = R. 

Recall that F is defined as F = $*0. Hence F* — 0*$. Since is two-sided inner, 
06* is identically equal to I r . Hence $ = OF*. 

It is well known that Tr > holds if and only if its outer spectral factor is 
invertible outer, c.f., [5J Proposition 10.2.1]. Assume Tr > 0. Then $ and 
are invertible with 3> _1 <E fT~v*» and _1 = 0*. This shows that F = $*0 is 

f/t A lit 

invertible in £„ xm , with inverse ($*0)" 1 = 9*($*) _1 = 6*($ -1 )*. Since 0* and 
(<j> -1 )* are both anti-analytic, so is F^ 1 . □ 

Proof of Theorem 10.21 We observed at the beginning of the present section that 
R admits an outer spectral factor and that is two-sided inner. The relations 
between the McMillan degrees of R, $, and F in (|0.11[) follow from Proposition 
Q The identity (|uT2"j) follows by replacing K in (JO]) by K = [ K F], noting that 
rank (TqTq - T R T* R ) = rank (TqT g - T K T* K - T F T*) < oo by Theorem and 
the identities T R T R = T K T* K + TpTp and H R H*~ = H K T£ + HpTp. □ 

In case (|3.2j) holds, the following proposition shows how the partial isometry M 
in (|0.6I) can be computed. 

Proposition 3.6. Let G G <RH™ xp and K e ^H^ xq such that ^372j holds. Let 
v = rank (TqTq — TkT^) < oo. Then v < S(K), the space Ho in (|0.5[) can be 
taken to be C v , and in that case the partial isometry M in (|0.5|) and (|0.6I) can be 
computed via M = MfM* with 

Mi = i / V(e iu )*V(e iw )dw, M* = ^- / V{e lul )*W{e luJ )duo, 
2?r J 2tt J q 

V(e lt ) = [ e lt A (e it ) G(e lt ) ] , W{e lt ) = [ A (e lt ) K(e lt ) ] , a.e. 
and Mi the Moore-Penrose pseudo-inverse of Mi. 

Proof. Recall from the introduction that dimHo = rank {T g Tq — T K T^) = v. 
Since v < oo, we can apply a linear transformation identifying Ho with C, and 
since T-L Q comes from the factorization of TqTq — T^T^, we can just as well apply 
this transformation and take Hq to be C y . The bound on v is a direct consequence 
of (EU). 

The formula for M Q follows by applying Lemma 14.21 with the given choice of V 
and W. Note that the identity (14.51) follows from (|0.4I) . The square summability of 
the Taylor coefficients of V and W follows from the boundedness of TcE p , TxE q 
and A (as defined in the introduction), as operators mapping into £+(C n ). Hence 
all conditions are satisfied, and Lemma 14.21 applies. □ 

We conclude this section with two examples. 

Example 3.7. According to Proposition 12.31 if H g Hq - H K H* K > 0, then the 
upper bound on the rank of TqTq — TkT^ — TpTp in item (i) can be improved 
to dim Ai , with Ai as defined in Lemma 12.11 This improvement can be arbitrarily 
large. Let I be a positive integer, take for G any rational function of McMillan 
degree I and take K = G. Clearly R = GG* - KK* = 0, thus $ = 0, which implies 
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A4 = {0}. Hence dim.M = 0; a solution with McMillan degree is obviously 
X(z) = 1, z £ C. On the other hand dim(Im Hq + Im Hk) = dim(Im Hq) — 
5(G) = I. Hence we have an improvement of /. 

Example 3.8. Let G and K are matrix polynomials whose values are matrices 
of size m x p, respectively m x q, say with degrees d±, respectively d 2 . Assume 
that the last coefficients of G and K, i.e, corresponding to z dl and z d2 , have full 
rank and that p,q> m. This implies that the last coefficients of G and K admit a 
right inverse. Note that Hq and Hk only have entries on the first d\, respectively 
d 2l anti-diagonals, starting in the left upper corner. Since the last coefficients 
of G and K admit a left inverse, it follows that 5(G) = ranking = m ■ d\ and 
5(K) = r&nkH K = m ■ d 2 . Now also assume that T G T G - T K T K > 0. Applying 
Theorem 10.11 and following the subsequent procedure we obtain that there exists 
a rational matrix solution A to (|0. 1 [) . The McMillan degree of X is bounded by 
5(G) + 5(K) — m(di + d 2 ). However, in this case the rank constraint in item 
(ii) of Theorem 10.11 gives a much sharper bound, namely dim(Im£f(3 + Im Hk) < 
m max{c?i, d 2 }, due to the specific structure of Hq and Hk- Note that this bound 
is in line with Q.3, (where the factor m does not appear, but should be there). 



In this appendix we prove two results of a general operator theoretical nature 
that are used in the paper. 

Lemma 4.1. Let V = V\ (B V 2 and W = Wi © W 2 be Hilbert space direct sums, 
and let X : V — > V and Y : W — > V be operators. Assume that X is selfadjoint and 
XV C Vi, and that Wi = Y~ [Vi], i.e., Wi is the inverse image of Vi under Y. 
Finally, let P-w t be the orthogonal projection of W onto Wi . Then 



Moreover, rank (YP^xY* — X) < dimVi- Assume YY* — X > and in addition 
that Y is injective and X > 0. Then rank (Y Pyy^*) = dim Wi, rank A < dim Wi 
and rank (YP Wl Y* — X) < dimWi. 

Proof. Using the decompositions V = V\ © V 2 and W = Wi © W 2 we represent A 
and Y as 2 x 2 operator matrices, as follows: 



(4.2) A = 



"Ai 0" 




"vr 




vr 


, Y = 


'Yi 


Y 2 




"Wi" 




"Vi" 







V 2_ 


V 2_ 





^3. 




w 2 _ 







Note that the zeros in the operator matrix for A follow from the fact that A is 
selfadjoint and XV C V\. The zero in the left lower corner of the operator matrix 
for Y is a consequence of Wi = y _1 [Vi]. Indeed, the latter equality implies that Y 
maps Wi into V%. The identity Wi = F _1 [Vi] also implies that Y 3 is one-to-one. 
To see this, assume Y 3 u = for some u £ W 2 . Then 7u 6 Vi. But the latter 
can only happen when u S y _1 [Vi] = Wi. Thus u G Wi fl W2, and hence u = 0. 
Therefore, I3 is one-to-one. 



4. Appendix 



(4.1) 




RATIONAL MATRIX SOLUTIONS TO THE LEECH EQUATION 



13 



Next, observe that the partitionings in (|4.2[) imply that 



(4.3) 



(4.4) 



YY* -X 







Y 2 
Y, 



Y{Y* -Xi 
Iw 2 



iVi u 

y 9 * Fo a 



Vi 

v 2 



FP W1 F 



FiYi* - Xi 




Vi 

V-2 



Now assume that the inequality in the right hand side of (|4.1|) holds. This 
implies that the operator matrix in the right hand side of (|4.4I) is positive. But 
then the same holds true for the operator defined by the second operator matrix 
in the right hand side of (|4.3|) . The equality (|4.3|l then shows that YY* - X is a 
positive operator, and the implication <= in (|4.ip is proved. 

To prove the reverse implication assume that YY* — X is a positive operator. 
Since I3 is one-to-one, the operator U from Vi © V2 to V\ ©W2 defined by the third 
operator matrix in the right hand side of (|4.3p has a dense range. Using (|4.4[) and 
the positivity of YY* — X, we see that 

YrYi* - X 1 



( 



C/w, C/w ) > for all v 6 V = Vi © V 2 . 



Iw 2 _ 

But the range of [/ is dense. Hence, by continuity, we get 
~Y X Y* -X l 



( 







y,y)>0 for all y e Vi © W 2 . 



It follows that YiYj* — Ai is positive, and by (|4.4I) the same holds true for the 
operator YPyj^Y* — X. This proves the implication =>■ in (|4.1|) . 

The decomposition (14. 4| shows clearly that rank (YPwjF* — X) < dim Vi. 

Note that if Y is injective, we have rank (YPyv^F*) — rank(Pwi) = dimWi. 
Assuming YY* — X > 0, we have YP Wl Y* > X. By Douglas' Factorization 
Lemma, = KPyy^Y* for some contraction K, and hence 

rank A = rank A ^ = rank {KP Wl Y*) < rank (P Wl F*) = rank (FP Wl F*). 

Thus rank X < dim Wi . A similar argument applied to FPyvi Y* > YP Wl Y * - X 
shows rank (YP m Y* — X) < rank (YP Wl Y * ) = dim Wi . □ 



Lemma 4.2. 

V(«) : C fe -> 



Consider two matrix functions V and W , analytic on D, with valued 
C p and W(z) : C u — > <C P . z G B, and Taylor expansions V(z) — 



J2^Lq z3 Vj an d W{z) — J^'jLo zJ Wj. Assume 



E%o VfVj < 00 and E~ WfWj < 



00. If 

(4.5) V(z)V(w)* = W(z)W(w)* for all z,w el, 

then there exists a partial isometry M : C — > C fc such that V(z)M 
z in B. Moreover, this partial isometry M is given by M 



W{z) for all 



M{NU with 



1 

M 1 = — 
2ir 



271 



V(e iu )*V(e iu )du 



OO 



(4.6) 



M* = 



2tt 



p27T OO 

/ V{e lu )*W{e lu} )dio = V V?W S 
J ° 3=0 



Here denotes the Moore-Penrose pseudo inverse of Mi 
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Proof. The assumption yields we can define operators fii and Q2 by 



:C k ^e 2 + (C p ) and n 2 



Wo 

w 2 



: C 



For each zeDwe write $ z for the point evaluation operator 



$ x = E%(I - zS;)- 1 :4(£7)^<C7, i.e., ^(x , x u x 2 , . . .) = ]T 
Note that V(z) = $ PtZ Sli and W(z) = $ p , z fl2, z £ D. Hence 

ff p ,»(nini - ^ns)^,™ - ^(2)^)' - - (*, w e d). 

Since n z6 DKer5"p iZ = {0}, it follows that fiifi* = f2 2 ^2 ■ By Douglas' factorization 
lemma there exists a unique partial isometry M : C' y — > C fe that satisfies OiM = O2 
and has Im fi^ as initial space and Im fi| as final space. Multiplying both sides with 
$ P ,z yields V(z)M — W(z), z £ D. Note that the Moore-Penrose pseudo inverse 
of 17 1 is given by 17^ = (f2^f2i) + f2 1 . Then 17+17i is the orthogonal projection on 
ImfiJ. Thus M = IJ+UjM = 17+17 2 = (nfni)+njn 2 - Note that Mi = fiffii and 
M* = njfi 2 - Hence M = M+M*. □ 
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